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1 Introduction 



The geometry of a manifold is given by a metric, which defines a notion of distance be- 
tween two points. Paths of shortest length connecting points are obtained as the critical 
curves of the functional variation of the integral defining arclength. Functional variation 
of this integral yields Euler-Lagrange equations which are a system of ordinary differen- 
tial equations of second order, whose solutions are the geodesies. Thus associated with 
geometry is a system of ODEs. This paper seeks to answer the inverse problem: when 
does a system of ODEs represent the paths of shortest length of a metric? That is, we 
wish to establish when ordinary differential equations exhibit an underlying geometry. We 
shall not be so ambitious as to attempt a solution on manifolds of arbitrary dimension 
and endowed with a general metric but shall restrict ourselves to the case of a pair of 
second-order ODEs on a (two-dimensional) surface and ask when the underlying geometry 
is flat, that is, a Euclidean space. We are concerned only with the solutions of the ODEs 
up to reparameterization since they serve merely to describe the paths of shortest length 
on the surface. Geodesies however, are not invariant with respect to general changes 
of parameter, so it shall be necessary to incorporate reparameterization in the precise 
definition of the problem. The formulation of the equivalence problem is the content of 
sections 1 and 2. In section 3, Cartan's method of equivalence is employed up to the 
level of the first normalization for generic ODEs. In section 4, we obtain the solution in 
the Euclidean case. The equivalence method is carried through in section 5, for generic 
pairs of second-order ODEs with the result that the symmetry of the system produces 24 
invariant functions. 

Before proceeding, it is intructive to consider the simpler problem without consid- 
erations of parameterization: when are the solutions to a pair of second-order ordinary 
differential equations 

d 2 y 1 /d 2 t = f(y,y,t) and d 2 y 2 jdh = g(y,y,t), (1) 

locally, the geodesies of some Euclidean metric on the plane? The geodesies of a Euclidean 
metric are the straight lines with respect to some coordinate system. Thus our problem 
may be formulated as follows: when does there exist a coordinate system Y = Y(y) such 
that the solutions y — y(t) of (1) correspond to straight lines Y = Y(y(t)) = at + b, a, b G 



R 2 ? We therefore seek to determine the existence of a transformation ^ : R 2 xR — + R 2 xR 
of the form 

such that \1/ transforms the equations 

d 2 Y 1 /d 2 t = and d 2 Y 2 /d 2 t = (2) 

into the equations (1). Any transformation of the form \I> above transforms (2) into 
equations of the form 

d 2 f/d 2 t + I ;;,(;/);>';/'' = 0. 

Consequently, (1) must necessarily be of this form. The terms r* fc (?/) define a connection 

V on the surface and we see that the solutions to (1) are locally the geodesies of a 
Euclidean space if and only if V is flat. The problem above requires that the solutions to 
(1) are already parameterized in such a fashion that only a change in the coordinates of 
the surface is sufficient to straighten them out into lines. In this paper, we are interested 
in whether the solutions to (1) may be reparameterized so as to be straight lines in some 
coordinate system. That is, do there exist coordinates Y — Y(y) and a reparameterization 
of time T = T(y,t) such that the solutions y = y(t) of (1) correspond to straight lines 

Y = Y{T) = aT + b, a,b G R 2 1 We therefore seek to determine the existence of a 
transformation $ : R 2 x R — > R 2 x R of the form 

*(y,t) = (Y(y),T(y,t)) 
such that $ transforms the equations 

d 2 Y l /d 2 T = and d 2 Y 2 /d 2 T = 

into the equations (1). 

Conceivably, other reparameterization criteria could be considered as well. For in- 
stance, one might investigate the more restricted transformation T = T(t) where time is 
reparameterized in a manner independent of the point on the surface or the more general 
T = T(y,y,t). Here, we shall content ourselves with spacetime reparameterizations only 
and defer the other cases to another time and place. 



It is convenient to place the problem in a more general setting; we consider the equiv- 
alence of two pairs of ordinary differential equations 

dV d 2 y 2 

= f(y, y, t) = g{y, y, t) (3) 

and 

d 2 Y x d 2 Y 2 



F(Y,Y,T) -— = G{Y,Y,T). (4) 



d 2 T v ' ' ' d 2 T 
under transformations of the form 

*(y,t) = (Y(y),T(y,t)). 

The case F = G = is solved in section 4. We obtain an e-structure on a 12-dimensional 
space with constant torsion. The solutions to a pair of ODEs belonging to this equivalence 
class have symmetries given by the group of fractal-linear transformations on the plane. 
In section 5 we make no restrictions on F, G, f and g and carry the equivalence through 
for the generic case. 

A similar problem was studied by S.S. Chern [3]. He has considered the geometry of 
a system of second order ODEs 

d 2 y i 

-jqr = f{y,y,t) i = l,...,n, 

under transformations of the form 

' Y = Y(y,t) 
T = t. 

Prior to Chern, the local behaviour of systems of second order ODEs has been studied by 
M.D.D. Kosambi [6] and by E. Cartan [2]. 



2 The Equivalence Problem Formulated 

The equations (3) may be represented by the Pfaffian system 



dy 1 


-p l dt 


= 


dy 2 


-p 2 dt 


= 


dp 1 


-fdt 


= 


dp 2 


— gdt 


= 



on U C R 5 . Similarly we represent the equations (4) by the Pfafhan system 



J 



dY 1 -P 1 dT = 

dY 2 -P 2 dT = 

dP l -FdT = 

dP 2 -GdT = 



on V C R b . 

Form the coframes 



where 





uj x = dy 1 




tt 1 = dY 1 




uj 2 = dy 1 + pdy 2 




fi 2 = dY l + PdY 2 


= < 


uj 3 = dy 2 + qdt and 


n = < 


Q 3 = dY 2 + QdT 




u) A = dp + hdt 




= dP + HdT 




u 5 = dq + gdt 

K 




VI 5 = dQ + GdT 


P 


= -P 1 (P 2 )~ 1 


p = 


-P\P 2 )- 1 


h 


= (p 2 )~ 2 (p 2 f - p 1 g) and * 


H = 


(P 2 )- 2 (P 2 F - P 1 ^) 


q 


= -p 2 


Q = 


-P 2 . 



Observe that / (resp. J) is spanned by uj 2 , ...,oj 5 (resp. f2 2 , ...,f2 5 ). 
Consider those transformations 

$ : U -> V 

such that 

' (1) $*(J) = I, and 

(*) 

(2) $(y,p,t) = (Y(y),P(y,p,t),T(y,t)). 
This is an overdetermined equivalence problem (cf. [5]). Nevertheless, we shall not follow 
the approach given in the above reference in that all the information contained in (*) may 
be encoded by an apprpriately chosen coframe and group, as follows. 

Let G be the subgroup of GL(5, R) whose elements are represented by 



/ 



a 


b 








c 














e / 







M 


N 



\ 



It is easily shown that a diffeomerphism $:[/—> V satisfies (*) if and only if 

for some 7 : U — > G. In order to avoid the awkward presence of the map 7, form the 
lifted coframe 77 (resp. H) on the G-bundle G x U (resp. G x V): 

77 := So;, (resp. if := S'fi), 

where 5 : G — > GL(5, R) is the natural injection. It can be shown that 

if and only if 

77, 

for some $:Gxf/-»GxF. i> is related to $ and 7 by <&(g,u) := (g^{u)~ l , for 
all E G x U. We arrive at the following formulation of the equivalence problem: 



The two systems of ordinary differential equations (3) and (4) are equivalent with 
respect to a diffeomorphism $ : U — > V of the form $(y,t) = (Y \y) , T (y , t)) if and only 
if there exists a diffeomorphism §(g,u) : G x U -»GxV satisfying 

$*(#) =77. 



With this characterization of equivalence we now proceed to reduce the group G. 

3 The First Normalization 

The structure equations, after making the obvious absorptions are 



a 


(3 





\ 




/ 





\ 





7 











Jrfrf + Aii] 1 !] 41 + A2fj 1 rj 5 








77 + 











e 










Bir] 1 r] 4: + B2 r q 1 rf 







n 


V 


) 




I 





/ 



The equations 

= d 2 n 2 n 2 
= d 2 rfrfrf 

give us the following infinitesimal group action on the torsion tensor. 



(5) 



dJ + (a -7 + 4>)J + A/i 2 = 
dA + (a - rfA + Av = > mod base. 

dB + {a- 4>)B + Bv-eA = 
Therefore the group action on the torsion tensor is 

J = {J Q -A Q N- l M 2 )a- 1 cf- 1 
A = a^cAoN- 1 
B = a-\fB + eA )N-\ 
where Jo, A , and B denote the tensors J, A and B, respectively, at the group identity. 
A parametric calculation shows that 

Jo = -ipq)-% A = (p~\ 0), B = (0, -(pgr 1 ). 

We may therefore normalize 

J = o, A = (0,1), and B = (1,0). 

Then 



t4 




= 


vl 




= 


e — 


v\ 


= 


a - 


1 + A 


= 


0" 


l-v\ + vl 


= 



mod base. 



Write 



Atf 
v\ = Brf 
e = vl + Ctf 
a = 7 — v\ + D(rf 

= 7 + v l - w 2 + 

Substituting these values into (5) results in 



B3 = At, D 3 = A 5 , D 4 = £? 5 , E 5 = C 4 + D 5 . 



The structure equations become 



drj 



7 — Vo 


Q 


o 


o 


o 


o 


'y 
i 


o 


o 


o 





v\ 7" 


- v\ - v 2 2 













iA 


v\ 


v\ 

















where 



T = 



A^rj 1 + B 5 r] 4 r] 1 + D 5 7] 5 7] 1 





Ci?? V + C 3 r] 3 r] 2 + C A r] 4 r] 2 + C b rfr] 2 + Eirfrf + E 2 r] 2 rf + E A tfif + (C 4 + D 5 )r] 5 7] 3 







Absorb the torsion as follows: 



Airfrf + A^rfrf + Bxrfrf + -E^ 5 ?] 

. v i + (£ 4 _ B 5 )rj 4 + C 4 rj 5 

^ ^ + (C 3 - £ 2 )?? 3 + C477 4 + C 5 ?? 5 

^ v 2 - A,?/ 3 - ^ - D 5V 5 



iA 



ii\ + (C 3 - E 2 ) v 5 



The new structure equations are 
/ 



drj 



1-vl 


P 














7 














vl 


7 + uj - 















v\ 


vl 
















V 

The equations 



rj + 



r} x rf + Arfrj 2 + Brfrf 







,1^4 



= d 2 rfifrf 

= (Prfrf^rf + dP-q^rfrfr) 4 

= d 2 rfrfr] 4 rf + d 2 rfrfr] 4 rf + d 2 rfr} 2 rfrf + d 2 rfr} 2 rfrf 



mod base. 



give the following infinitesimal group action on the torsion tensor: 

dA + A(7 -v{) + Bvl + 2(j,{ = 

dB + (7 - u*).B + 2v\D + 2/4 -2/^ = 

dC + (2 7 + v\ - 3uf)C + n\D = 

dD + (7 + uj - 2u|)D = 

4 Geodesies of Flat, Symmetric Connections 

In this section the equivalence problem is carried through for the case F = G = 0. This 
will lead to an e-structure with constant torsion on a 12-dimensional space. The only 
invariants are therefore constant invariants and thus we obtain a complete solution to the 
problem of equivalence. 



4.1 The Second Normalization 

A parametric calculation will show that at the identity A = B = C = D = 0, hence 
D = C = 0. This leaves the following two equations: 



dA + A(7 - v\) + Bv\ + 2fi\ = 
dB + (7 - + 2/4 -2/^ = ^ 

Normalize A = -B = 0. We then have 



mod base. 



/4 







mod base. 



This produces new torsion by 

This gives the following structure equations: 



1A 



drj = 





P 














7 














v\ 7" 


- v\ - v\ 
















vl 


v\ 

















+ T, 



where 



T = 







A 1 i] 1 i] 2 + A^rfrf + A A rfrf + A 5 i] 5 i] 2 + B^i] 3 + B 2 i] 2 i] 3 + B 4 i] 4 i] 3 + B b rfrf 





Absorb 



We then obtain 
( 



dr) 



V 

The equations 



vl — > v I - A 4 r] 2 - B 4 rf 

v\ — > f 2 — A 5 r] 2 — _B 5 7y 3 

At? -> //? - A 3 r] 2 + S 2 ?7 2 

7 -> 7 + A^ 2 - .Bgr] 2 

/3 -> (3 + A^-B^ 1 . 



\ 



7-^2 

















7 
















7 + uj - u| 























Mi 











1] + 




Arfrf + Brfrf 




= (PrflfrfTf 

= d^rf^rf + AV^V 



give 



cL4 + (27 - uj - u|)i4 + 5^ = 
dB + (2~f - 2v%)B = 

Now at the identity A = -Bo = 0. Thus A = B — 0, and hence 



mod base. 



( 



drj 



V 



7-^2 

















7 
















7 + vl - uf 














Mi 


vl 


vl 
















\ 




Only constant torsion remains so the system must be prolonged. 



4.2 Prolongation 

Now dim = 2. The first prolongation corresponds to the following arbitrariness in 
the tableau: 

9 1 := v\ 

9 2 := v\ + at] 3 

9 3 := v\ + arf 

9 4 := 1 + 2arf 

9 5 ■= fi + ar] 1 + br] 2 

9 6 := ix\ + arf 

where a, b e -R are arbitrary. 

The structure equations for r\ may be written, 



9 4 - e 3 


9 5 











\ 




( 








9 4 

















1 ^ 
7] 7] 





e 2 


e A + e 1 - e 3 










7] + 




r] l r] 4 










e 1 


e 1 






















e 3 


J 




K 






Taking d 2 of the above structure equations we obtain the following: 

o = -de 3 ?] 1 + c$y + d0V + 3 0V - 

= d9 4 rj 2 + 9 5 rj 2 rj 5 + 9 6 rj 1 rj 2 

= dfly + ^2 _ ^3 + ^3 _ 01^2 + £3^2 + £5^4 + 

= d9 1 r ] 4 + c/6> 2 77 5 + dO^rf - 9 2 9 3 V 5 + 9 6 9 3 V 3 + 9 4 9 6 V 3 - 9 6 r ] 1 r ] 4 - 9 1 9 2 r ] 5 

o = d9 3 v 5 + de 6 r] 2 - 9 3 9 6 r] 2 + e A e 6 v 2 -e 6 ^ 5 . 

It follows from a somewhat lengthy but straightforward calculation that the structure 
equations for 9 are given by 

d9 1 = Airfrf 2 + A 2 r] 3 r] 2 + A377V + A^fif + A577V + 9 b rf> + 6»V 

^2 = $ y + Air] A v 2 + £^2 + ^4^3 + ^2 _ £3^2 + £5^4 

^3 = ^ r j^ + Crj 5 T] 2 + 29 5 rj 5 + 9 6 rj 1 

d9 4 = 2$y + D x ifi] 2 + z^y + (c - + #y - #y 

^5 = $ y + $2^2 + + D 2 r ] 4 r] 1 - A x ifi] x - 9 3 9 5 

d 06 = $ y + prftf + Aa?7 Y + 0306 _ QAQ^ 



where $ 1 and $ 2 are 1-forms in the new tableau. The equations 

= d 2 6 1 

= d 2 e 2 v 3 

= d 2 6 2 rf + d 2 6 z rf 

o = 2d 2 e 3 -d 2 e 4 

= d 2 Q % rf 

give the following infinitesimal action on the torsion tensor 



dA 1 - $ 2 

dA 2 

dA 3 

dA A 

dA 5 

dB 












> mod base. 



dC-2® 2 = 

dD 1 = 

dD 2 = 

dF = 



At the group identity, 



A 2 = A 3 = A A = A 5 = B = D l = D 2 = F = 0, 



therefore these terms are identically zero. Thus, 



Consequently, 



We also have, 



Therefore, 



(dA 1 + A 1 (6 S + 4 ) - ^j^V + (2A 1 - C)e 2 7] 5 7] 2 7] 6 = 



C = 1A X . 



dA 1 + Ai(0 3 + 4 ) - $ 2 = mod ?7 2 , 77 5 



cL4i + Ai(0 d + 4 ) - $ 2 + V = 



for some function /. We obtain the following structure equations 

dO 1 = Arfrj 2 + 0V + 0V 

d9 2 = $1^3 + ^4^2 + 0102 _ q3q2 + 9 5 v 4 

d9 3 = $ y + 2Arfrf + 20V + #V 
d0 4 = 2$V + Arfrf + 0V - 6 ?? 1 

^5 = $y + $2^2 _ ^5^1 _ 3 5 

where A has been written for Ai. 

4.3 The Third Normalization 

We have the infinitesimal group action on the torsion given by 

<L4 1 + A 1 (0 3 + 4 )-$ 2 = O mod base. 
Normalize A = 0. Then $ 2 = Z?? 2 and hence 

de 1 = 0V + ^V 

rf0 2 = $1^3 + 0l#2 _ 9 3 9 2 + ^4 
rf03 = $ 1 r ? 2 + 20V + ^V 

rf0 4 = 2$v + ^V-^V 

d05 = $y - 03^5 

d06 = $1^5 + 0306 _ 0406 

We have constant torsion and an e-structure. Thus (77 1 , rf , 1 , 6 , $) is an invariant 
coframe (here we have written $ = .) The equations 

(d$ - $0 4 - 5 6 )^ = 

for % — 1, 2, 3 and 5 give 

d$ = $0 4 + 0^6. 



Thus the structure equations are 



dr] 1 
d V 2 
drf 
drf 
drf 
d6 l 

de 2 

d6 3 

de 4 
de 5 
de 6 



(e 4 - e 3 )^ 1 + #y 
#y + rfrf 

e 2 v 2 + {e 4 + e l - e 3 )rf + rfrf 
0V + #y + e 2 rf 
e 6 v 2 + e 3 rf 
#y + #y 

$r? 3 + (6 l - 6 3 )6 2 + e 5 v 4 
<&rf + 29 b rf + 6 & rf 
2<$>rf + #V - Q % rf 

<5>ri 5 + (# 3 - # 4 )# 6 

<&e 4 + # 5 # 6 . 



We have shown the following: 



Theorem 1 The solutions y(t) of 



(6) 



d y /d t = f(y,y,t) and d y jd t = g(y,y,t), 



(7) 



are the geodesies of a Euclidean space with respect to a transformation of the form (Y, T) = 
(Y(y),T(y,t)) if and only if equations (7) yield the structure equations (6). 

The structure equations (6) are the Mauer-Cartan equations for a Lie-group: the 
group of fractal-linear transformations T on the plane. Recall that the fractal-linear 
transformations on the plane are those transformations A of the form 

A = (y,i) = (y{y),i(y,t)) 



where 



y\y) = 
y\y) = 
t(y,t) = 



bl + bly 1 + h\y 2 
a + a x y l + a 2 y 2 
b 2 + bjy 1 + b\y 2 



a + a x y 4 + a 2 y 2 



t + c + ciy 1 



c 2 y 



a + aty 1 + a 2 y 2 



where a*, a*-, 6*- £ R are constants. Let Q denote the subgroup of GL(A,R) consisting 
of those invertible matrices M whose first column is '(1,0,0,0). We may identify the 
transformation A with the element <p{A) in Q by 





c 


Cl 


c 2 


\ 





a 


a\ 


a 2 







b 1 


b\ 


b 1 







b 2 


bi 


b 2 


J 



The map : T — > Q is a group isomorphism and so we may identify the fractal-linear 
transformations with Q. 



5 The Generic Case 

In this section we continue the equivalence problem from the end of the first normalization 
(section 3). 



5.1 The Second Normalization 



Recall that after the first normalization we have the structure equations 





f 7 


-v\ 


p 


















7 











drj = 







v\ 


7 + v\ - v\ 















iA 


A. 


v\ 


v\ 




I 













v\ 



) 



1 ^ 

r] < rj 



i] 1 !] 4 + Arfrj 2 + Br] l r] 



4^3 



with infinitesimal group action on the torsion tensor given by 

dA + A(7 - v\) + Bu£ + 2fj\ =0 

cLB + (7 - vl)B + 2v\D + 2/4 - 2/4 = 

dC + (2 7 + uj - 3v 2 )C + \i\D = 

dD + (7 + v{ - 2vl)D = 

In the generic case, We thus normalize 







^ Cr] 1 !] 3 + Drfrf J 



mod base. 



A = 0, 5 = 0, C = 0, D = 1. 



This gives us 



iA 



o 
o 



7 = 2vj-v\ 



mod base. 



Let 



lA 

iA 
tA 



7 = 



It follows that 

The new structure equations are: 
/ 



vl + Ctf 

1v\ - v I + Drf 

B d = D 3 . 



dr) 



vl - v\ 


P 















/ 










2vl - vl 

















i ^ 

7] 7] 







vl 


vl 










rj + 




1 4 
7] 7] 













vl 


vl 
















vl 








vl 


) 




I 


1 4 
7] 7] 


) 



+ T, 



where 



Btffv) 1 + D 4 7] 4 7] 1 + D^TfT] 1 + D 2 77 2 ?7 



5„1 



-,2^1 



T = 



D^T] 2 + B^T] 2 + D 4 7] A 7] 2 + D 5 7] 5 7] 2 

D 1 7 ] 1 7] 3 + I^V + D^V + D 5 ??y 



Ai?] 1 ^ 2 + A 3 7] 3 7] 2 + A 4 7] A 7] 2 + A 5 7] 5 7] 2 + Bi7] l 7f + B 2 7] 2 7] 3 + B^T] 3 + B 5 7] 5 7} 3 

C^T] 2 + C 3 7?V + C477V + C 5 r?V 



Absorb the torsion: 







uj - 5 4 ?f - D 4 ?? 4 - A 4 r? 2 

u£ - A 5 7] 2 + C 3 7] 3 

(3 + A^ 1 + C377 1 
t; 2 + (C 3 + D 2 )r ? 2 . 



The structure equations are then 



drj 



Vn — v} 


ft 
p 


n 


n 


n 


o 


o v 2 — v } 


o 


o 


o 





vl 


vl 

















v\ 


vl 





v\ 











V 

5.2 The Third Normalization 

The equations 



r]+ 



Arftf 



77V + Bifif + Arfrf 
i] 1 i] A + Bi] 1 i] 3 + + Arfrf 

Drj 1 !] 2 + Erfr} 2 + F?/?? 3 + Grfrf 



= 


(Prfrfrf + (Prfrfrf -f 


d 2 7] s r] 1 r] 2 


= 


d 2 rfrfrf + cPrfrfrf 




= 


dPrfrfrf 




= 


d 2 ifr) 2r if — d 2 rfri 4 rf 




= 


d 2 rfrfrf - d 2 rfrfrf 




= 


-d 2 rfr] A rf - d 2 rf , rfrf 


+ d 2 ifr\ 2 if + d 2 r] 5 r] 3 r] 



(8) 



give us the following infinitesimal group action on the torsion tensor: 



dA + Av\ + Cvl + 2(3 
dB + B{v 2 -v\) 
dC + Cv\ 

dD + 3D(v% -v\) + Fv\ - Hv\ 

dE + E(3v% - 2v\) — F(3 — Gv\ 

dF + 2F(v 2 -v{) 

dG + G(2v 2 -v\) 

dH + 2H(v 2 - vl) - Bvl 

dl + Iv 2 - Cvl 

dJ + J(2v\ - vl) + B(3 + Iv\ 
















mod base. 



At this point, the group is the subgroup of GL(5, R) whose elenents are of the form 



ca 1 


e 











c^- 1 











a^bc c 














a 


b 





or x bc 





c 



The group action on the torsion is then given by 

A = AqC 1 - Coa^bc' 1 - 2ac~ 2 e 
B = Boac- 1 
C = Coa- 1 

D = D a 3 c- 3 + (l/2)B ab 2 c- 3 + (H -F )a 2 bc- 3 

E = E a 2 c~ 3 + F a 3 c" 4 e + G afec- 3 

F = F a 2 c- 2 

G = G ac~ 2 

H = H a 2 c~ 2 + B abc- 2 

I = / oC -i + Coa^bc- 1 

J = J ac~ 2 -B a 2 c- 3 e-{l/2)C a- 1 b 2 c~ 2 -I bc- 2 . 
We may then normalize 

A = 0, B = l, C=l, 1 = 0. 

It follows that 

v\ = v\ = v\ = P = mod base. 

We write 

v\ = Atf 

V 2 = BiT] 1 

v\ = C lV l 
P = D iV \ 

Substituting these values back into the equations (8) we obtain 

A 5 = 2D A + C 4 and .64 = ^1 + A 4 . 

The group has been reduced to the identity and so we have an e-structure. The structure 
equations are 



dr/ 1 = (£>! + A 2 - B 2 )r ] 1 r] 2 + (A 3 - B^rfrf - A^r] 4 + (A 5 - B^r] 5 - D 3 r] 2 r] 3 
+ (1/2)(C 4 -A 5 )77V-^V 

drj 2 = (1 + 25! - A 1 )r ] 1 r ] 2 + 77V + (A 3 - 2B 3 )r ] 2 r ] 3 - (A, + B 4 ) V 2 V 4 + (A 5 - 2B 5 

dr/ 3 = drfrf + (1 + B^rfrf + rfif + (S 2 - C 3 )r\ 2 if - C477V - C^rfrf 
- (1 + B 4 )»7V - £577 V 

dr] 4 = I 1 i] 1 i] 2 + I 3 i] 1 i] 3 + Axi] 1 '!] 4 + C 1 i] 1 i] 5 — I 2 i] 2 f] 3 + A 2 rfr] 4 + C 2 rfrf + A 3 i] 3 i] 4 
+ (C 3 - / 4 )^7 3 r7 5 + (C 4 - A, W 

cfy 5 = (Ci + / 5 )?/V + 7?V + Bi^V - C 4 r?y + (#2 - / 6 - C^V + B 3 r] 3 r] 5 
+ B A r] 4 rf - C 3 r] 2 rj 3 . 



We obtain 24 local invariants. 
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